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Protein lysine malonylation, a newly identified protein 
post-translational modification (PTM), has been proved to 
be evolutionarily conserved and is present in both eu- 
karyotic and prokaryotic cells. However, its potential 
roles associated with human diseases remain largely 
unknown. In the present study, we observed an elevated 
lysine malonylation in a screening of seven lysine acy- 
lations in liver tissues of db/db mice, which is a typical 
model of type 2 diabetes. We also detected an elevated 
lysine malonylation in ob/ob mice, which is another 
model of type 2 diabetes. We then performed affinity 
enrichment coupled with proteomic analysis on liver 
tissues of both wild-type (wt) and db/db mice and iden- 
tified a total of 573 malonylated lysine sites from 268 
proteins. There were more malonylated lysine sites and 
proteins in db/db than in wt mice. Five proteins with 
elevated malonylation were verified by immunoprecipi- 
tation coupled with Western blot analysis. Bioinformatic 
analysis of the proteomic results revealed the enrich- 
ment of malonylated proteins in metabolic pathways, 
especially those involved in glucose and fatty acid 
metabolism. In addition, the biological role of lysine 
malonylation was validated in an enzyme of the glycol- 
ysis pathway. Together, our findings support a potential 
role of protein lysine malonylation in type 2 diabetes 
with possible implications for its therapy in the 
future. Molecular & Cellular Proteomics 14: 10.1074/ 
mcep.M114.041947, 227-236, 2015. 
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Post-translational modifications (PTMs)' have been recog- 
nized as a common feature of proteins (1-3). More than 300 
types of PTMs have been identified according to the Swiss-Prot 
database (4, 5). Most of them use small molecular compounds 
as group donors. For example, adenosine-triphosphate (ATP) is 
used in phosphorylation, S-adenosylmethionine (SAM) in meth- 
ylation, and acetyl-CoA in acetylation. Lysine acylations includ- 
ing malonylation (6), succinylation (7), butyrylation (8), propiony- 
lation (9), and crotonylation (10) represent a group of PTMs that 
use intermediates of energy metabolism like malonyl-CoA, suc- 
cinyl-CoA, butyryl-CoA, propionyl-CoA, and crotonyl-CoA as 
group donors. Among the lysine acylations, lysine malonylation 
was first identified in Escherichia coli (E. coli) and HeLa cells 
using a specific anti-Kmal (anti-malonyllysine) antibody (6). It 
was found in three proteins in E. coli and 17 proteins in HeLa 
cells. Using a novel chemical fluorescent probe, another group 
identified more than 300 malonylated protein candidates in 
HeLa cells (11). Despite the rapid progress in detection tech- 
nologies and tools, functional studies of lysine malonylation and 
its role in human diseases have been lagging behind. 

Type 2 diabetes is characterized by hyperglycemia and 
production of glycated proteins. For example, glycated he- 
moglobin A1c (HbA1c) has been clinically used as diagnostic 
criteria for diabetes. In addition to glycation, the role of other 
types of PTMs in type 2 diabetes remains to be revealed. In 
fact, elevated malonyl-CoA levels have been found in type 2 
diabetic patients (12), and prediabetic rats (13). And hepatic 
overexpression of malonyl-CoA decarboxylase (MCD) de- 
creased malonyl-CoA and reversed insulin resistance (14). 
Given the use of malonyl-CoA as malonyl donor in lysine 
malonylation, lysine malonylation is therefore anticipated to 
be of functional significance in the pathogenesis of type 2 
diabetes. 


1 The abbreviations used are: PTM, post-translational modifications; 
10-FTHFDH, 10-Formyltetrahydrofolate Dehydrogenase; ACAT1, 
Acetoacetyl-CoA thiolase; ALDOB, fructose bisphosphate aldolase 
B; FBP1, fructose-1,6-bisphosphatase 1; G6PI, glucose-6-phos- 
phate isomerase; GSTs, glutathione S-transferases; HbA1c, hae- 
moglobin A1c; HMGCS2, Hydroxymethylglutaryl-CoA synthase, 
mitochondrial; LDHA, lactate dehydrogenase A; MCD, malonyl- 
CoA decarboxylase; PBE, Peroxisomal bifunctional enzyme; SAM, 
S-adenosylmethionine. 
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In the present study, we observed elevated lysine malony- 
lation in liver tissues of db/db mice after unbiased screening 
seven types of lysine acylations. We then detected elevated 
levels of lysine malonylation in liver tissues of more db/db and 
ob/ob mice. Using an immunoaffinity based proteomic 
method, we identified a total of 573 malonylated lysine sites 
from 268 proteins in liver tissues of wt and db/db mice. 
Elevation of lysine malonylation in five proteins was confirmed 
by immunoprecipitation coupled with Western blot analysis. 
Functional analysis of the malonylated proteins showed an 
apparent enrichment in metabolic pathways, especially those 
involved in the glucose and fatty acid metabolism. Our study 
indicates the putative association between protein lysine ma- 
lonylation and type 2 diabetes. 


EXPERIMENTAL PROCEDURES 


Reagents—Anti-malonyllysine (PTM-901), anti-acetyllysine (PTM- 
105), anti-1,2-dimethyllysine (PTM-602), anti-succinyllysine (PTM-401), 
anti-butyryllysine (PTM-301), anti-propionyllysine (PTM-201), anti-cro- 
tonyllysine (PTM-502) antibodies, and anti-malonyllysine conjugated 
agarose beads were purchased from PTM Biolabs (Chicago, IL). G6PI 
(ab86950), LDHA (ab47010), FBP1 (ab109732), and 10-FTHFDH/ 
ALDH1L1 (ab56777) antibodies were purchased from Abcam (Cam- 
bridge, MA). Glutathione (GSH) agarose beads were purchased from GE 
Healthcare Life Sciences (Uppsala, Sweden). Sequencing-grade trypsin 
was purchased from Promega (Madison, WI), and C18 ZipTips were 
purchased from Millipore (Billerica, MA). A list of other chemicals used 
for LC-MS/MS is available in previously published study (15). 

Animals and Samples—C57 BLKS db/db mice (n = 5) and wt mice 
(n = 5) were provided by the Animal Center of Peking University. C57 
BL/6J-ob/ob mice (n = 4) and wt mice (n = 4) were provided by the 
Institute of Laboratory Animal Science, CAMS (Chinese Academy of 
Medical Science). Mouse studies were approved by the Animal Ex- 
perimentation Ethics Committee (Institute of Biophysics, CAS) and 
the National Health and Medical Research Council of China Guide- 
lines on Animal Experimentation. 

Tissue Lysate Preparation— Mice were sacrificed by cervical dislo- 
cation. Tissues were homogenized in 1 ml lysis buffer (50 mm Tris-Cl, 
pH 7.4, 150 mm NaCl, 0.5 mm EDTA, 1 mm DTT, 1% Triton X-100, 
0.5% sodium deoxycholate, and 0.1% SDS). The homogenized tis- 
sues were centrifuged at 12,000 rpm and the supernatant containing 
the tissue lysates was kept for later steps. 

Western blotting and Coomassie Blue Staining—Western blot and 
Coomassie blue staining analysis were carried out as described pre- 
viously (16). 

In Vitro Malonylation Reaction and ALDOB Activity Assay—lIn vitro 
malonylation was performed as previously described (17). Briefly, 
FLAG-tagged ALDOB was overexpressed in 293T cells and purified 
using ANTI-FLAG M2 Affinity Gel (A2220, Sigma-Aldrich, St. Louis, 
MO) as per instruction by the manufacturer. Purified ALDOB proteins 
were incubated for 12 h at 37 °C in buffer containing 1 mm malonyl- 
CoA, 50 mm Tris-Cl, pH 8.0, and 150 mm NaCl. 

Enzymatic activity of ALDOB was determined spectrophotometri- 
cally by measuring the decrease of NADH at 340 nm in a coupled 
assay with triosephosphate isomerase (TPI) and glycerophosphate 
dehydrogenase (GDH) (18). Assays were performed in triplicates at 
30 °C in buffer A (100 mm Tris-Cl, pH 7.5, and 1 mm EDTA) following 
addition of fructose 1,6-bisphosphate. 

In-gel and In-solution Tryptic Digestion—For in-gel digestion, pro- 
tein bands around 95 KDa, 55 KDa, 34 KDa, and 26 KDa were excised 
from the SDS-PAGE gels, cut into small pieces and destained using 


50 mm triethyl ammonium bicarbonate (TEAB) in 40% acetonitrile. 
100% acetonitrile was used to dehydrate the gel. 5 mm Tris (2-carboxy- 
ethyl) phosphine (TCEP) was used to reduce the proteins. 10 mm methyl 
methanethiosulfonate (MMTS) was used to alkylate the proteins. Diges- 
tion was carried out using sequencing-grade modified trypsin (Pro- 
mega) with a substrate ratio of 1:50 (w/w) overnight at 37 °C. Peptides 
were extracted with 60% acetonitrile and 5% formic acid (FA). Resulting 
peptides were dried and prepared for immunoprecipitation. For in- 
solution digestion, liver tissue lysates were precipitated with acetone 
and resolved in 8 m urea with 50 mm Tris-HCl, pH 8.0. The samples were 
reduced and alkylated as done for in-gel digestion. 

Affinity Enrichment of Malonylated Peptides—The tryptic peptides 
obtained from in-gel digestion or in-solution digestion were resolubi- 
lized in NETN buffer (60 mm Tris-Cl, pH 8.0, 0.5% Nonidet P-40, 100 mm 
NaCl, and 1 mm EDTA). Anti-malonyllysine antibody-conjugated aga- 
rose beads were added and incubated at 4 °C for 12 h with gentle 
shaking. The beads were washed two times with NETN buffer, once 
with ETN buffer (50 mm Tris-Cl, pH 8.0, 100 mm NaCl and 1 mm EDTA), 
and once with purified water. The bound peptides were eluted by 0.1% 
trifluoroacetic acid (TFA) and dried using the SpeedVac. The dried 
peptides were re-dissolved in 0.1% FA and desalted using C18 ZipTips. 

Protein Identification and Quantification by LC-MS/MS—LC- 
MS/MS analyses were performed on a LTQ-Orbitrap mass spectrom- 
eter (Thermo, San Jose, CA) equipped with an Eksigent nanoLC 
system (Eksigent Technologies, LLC, Dublin, CA). The peptide mix- 
tures were loaded onto 20-mm ReproSil-Pur C18-AQ (Dr. Maisch 
GmbH, Ammerbuch) trapping columns (packed in-house, i.d., 150 
pm; resin, 5 wm) and eluted into 150 mm ReproSil-Pur C18-AQ (Dr. 
Maisch GmbH, Ammerbuch) analytical columns (packed in-house, 
i.d., 75 um; resin, 3 um) for MS analysis. Elution was achieved with a 
0-35% gradient buffer B (Buffer A, 0.1% formic acid and 5% aceto- 
nitrile; Buffer B, 0.1% formic acid and 95% acetonitrile) over 90 min 
using a data-dependent method. All MS/MS spectra were collected 
using normalized collision energy (setting, 35%), an isolation window 
of 3 m/z, and one micro-scan. The acquired MS/MS spectra were 
searched using the Mascot (v2.3.2) search engine in the Proteome 
Discoverer™ Software (v1.3, Thermo Fisher Scientific, Bremen, Ger- 
many) against the mouse Uniprot database (version 2.3, 50,807 se- 
quences) using a target-decoy database searching strategy (19). The 
following search criteria were employed: three missed cleavages 
were allowed; mass tolerance for parent ions was set to 20 ppm; 
mass tolerance for fragment ions was set to 0.8 Da; minimum peptide 
length was set at 6; Cys (+45.9877 Da, methanethiosulfonate, MMTS) 
was set as a fixed modification; Met (+ 15.9949 Da, oxidation) and Lys 
(+86.0004 Da, malonylation) were considered as variable modifica- 
tions. The global false discovery rate (FDR) for both peptides and 
proteins was set to 0.01. The relative quantitation of proteins between 
two samples (wt versus db/db mice) was achieved by extracted ion 
chromatograms (XICs) of peptides with a Mascot score of =25 that 
ranked first. The resulting peptides were extracted using the Precur- 
sor lon Area Detector feature in the Proteome Discoverer program 
(version 1.3), using a mass tolerance of 3 ppm. Only proteins with a 
minimum of two quantifiable peptides were included in our final data 
set. The peptides did not match any proteins present in the Uniprot 
database, and peptides with no XICs in both wt and db/db mice were 
excluded from the final data set. 

Bioinformatics Analysis of Malonylated Proteins—The malonylated 
proteins were classified based on the PANTHER (Protein Analysis 
through Evolutionary Relationships) system. The DAVID 6.7 was em- 
ployed for enrichment analysis of Gene Ontology (GO) terms and 
Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The 
default mouse proteome was used as the background list. The sig- 
nificance of the obtained enrichments was evaluated by using the 
Benjamini-Hochberg corrected p value. The identified proteins with 
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lysine malonylation were searched against the STRING database (V 
9.1) for protein-protein interactions. A network of protein-protein 
interactions was generated, which was then visualized using the 
Cytoscape program (v3.0.2). This interaction network was further 
analyzed for densely connected regions using the program “Molec- 
ular Complex Detection” (MCODE). 


RESULTS 


Protein Lysine Malonylation Elevate in Liver Tissues of 
db/db and ob/ob Mice—The db/db mice, which has a ho- 
mozygous point mutation in the db gene encoding for the 
leptin receptor, is an classical type 2 diabetic mouse model 
(20, 21). To investigate the role of protein lysine acylations in 
type 2 diabetes, seven types of lysine acylations were inves- 
tigated by Western blot in liver tissue lysates of wt and db/db 
mice. A marked elevation of lysine malonylation was observed 
in db/db mice relative to wild type littermates (Fig. 1A), while 
little elevation of lysine acetylation, 1, 2-methylation, succiny- 
lation, butyrylation, propionylation, and crotonylation was de- 
tected (Fig. 1B-1G). Equal loading was verified by B-actin 
Western blot and Coomassie blue staining (Fig. 1H, /). Fur- 
thermore, we detected lysine malonylation in another four wt 


wt ob wt ob 


Fic. 2. Detection of lysine malonyla- 
tion in liver of four db/db and ob/ob 
mice. Western blot analysis of liver tis- 
sue lysates from four db/db (db) A, and 
ob/ob (ob) B, mice and their counterpart 
wt mice by anti-malonyllysine antibody. 
Equal loading was verified using B-actin 
Western blot and Coomassie blue 
staining. 


B-actin malonyl 


coomassie 


and db/db mice and observed a significant increase in lysine 
malonylation in liver tissues of all db/db mice (Fig. 2A). To 
investigate whether lysine malonylation levels are elevated 
beyond wt levels in other type 2 diabetic animal models as 
well, we tested lysine malonylation levels in ob/ob mice. Con- 
sistent with the observations in db/db mice, an elevated levels 
of lysine malonylation was detected in liver tissues of all ob/ob 
mice relative to wt mice (Fig. 2B). Together, these findings 
demonstrated that lysine malonylation was elevated in liver 
tissues of db/db and ob/ob type 2 diabetic animal models. 
Identification of Malonylated Proteins in Gels—To validate 
the increase of protein lysine malonylation and identify the 
malonylated lysine sites and proteins, liver lysates from three 
wt or db/db mice were combined and proteins were sepa- 
rated by SDS-PAGE. Four gel bands with the highest levels of 
lysine malonylation were cut and trypsin-digested (Fig. 3A). 
For analysis of each band, a total of 50 gels were cut and 
combined. The digested peptides were immunoprecipitated 
with anti-malonyllysine antibody-conjugated agarose beads 
and then loaded onto LC-MS/MS. One representative MS/MS 
spectrum of the identified malonylated peptides is shown in 
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Fic. 3. Identification of protein lysine malonylation by in-gel digestion and LC-MS/MS. A, Diagram of in-gel digestion coupled with 
LC-MS/MS analysis. Four pairs of bands in Coomassie blue stained gel, which showed great increased levels of lysine malonylation by Western 
blot, were cut and trypsin-digested in-gel. Digested peptides were extracted and immunoprecipitated using anti-malonyllysine antibody- 
conjugated agarose beads. Immunoprecipitated peptides were identified by LC-MS/MS. x50: a total of 50 gels were cut for each bands. B, 
A representative MS/MS spectra of malonated peptides. The neutral loss of 44KDa is unique signature of malonylation and designated as 
“—44.” C, Number of malonylated proteins and peptides detected for each band. 


Fig. 3B. There was a neutral loss of carbon dioxide (CO., 44 
Da), which was used for confirming the presence of lysine 
malonylation. The LC-MS/MS results showed that the number 
of malonylated proteins and peptides identified for each of the 
four bands was higher in db/db mice relative to wt (Fig. 3C). 
Detailed information on all malonylated proteins and peptides 
identified is summarized in supplemental Table S1. Together, 
our results indicated that proteins from the livers of db/db mice 


were more intensively malonylated than those in wt, which was 
consistent with the results obtained from Western blot analysis. 

Identification of Malonylated Proteins in Solution—Having 
successfully tested the ability of our LC-MS/MS method to 
detect malonylation, we used in-solution tryptic digestion to 
widen the scope of our analysis and identify lysine malonyla- 
tion in a larger number of peptides (Fig. 4A). Liver lysates were 
combined from either three wt or three db/db mice, then 
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Fic. 4. Proteomic identification and 
bioinformatics analysis of protein ly- 
sine malonylation. A, Diagram of in-so- 
lution digestion coupled with LC-MS/MS 
analysis. Liver tissues from wt and db/db 
(db) mice were lysed. Proteins were 
extracted, trysin-digested, immunopre- 
cipitated using anti-malonyllysine anti- 
body-conjugated agarose beads, and 
detected by LC-MS/MS. x3: livers tis- 
sues from three mice were combined. B, 
Venn diagram showing overlapped ma- 
lonylated peptides between wt and db. 
C, Venn diagram showing overlapped 
malonylated proteins between wt and 
db. D, Pie diagram showing biological 
process enrichment of total malonylated 
proteins. Analysis was done with the 
web tool “PANTHER.” E, Network dia- 
gram showing interactions between 
malonylated proteins extracted from 
STRING database. MCODE method was 
used to calculate clusters. Five clusters 
were identified. 
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Fatty acid metabolism 


trypsin-digested and immunoprecipitated using anti-malonyl- 
lysine-conjugated agarose beads. After immunoprecipitation, 
LC-MS/MS, and a search against the mouse Uniprot data- 
base, we identified a total of 268 Uniprot-reviewed proteins, 
containing 573 malonylated lysine sites (supplemental Table 
S2). There were 246 lysine malonylation sites that overlapped 
between wt and db/db mice (Fig. 48); and 169 malonylated 
proteins that overlapped between the two groups (Fig. 4C). 
We extracted ion chromatograms and calculated the area 
under the peak associated with malonylated peptides, which 
indirectly reflects the quantity of each malonylated peptide. 
For the 246 overlapping peptides, the average under-the- 
curve area of extracted ion chromatograms for each malony- 
lated peptides in db/db was larger than that in the wt mice 
(supplemental Fig. S1). Together, these results suggested 
greater malonylation in db/db livers relative to wt, and were 
consistent with the results obtained by Western blot and 
in-gel identification. 


Glucose metabolism 


We compared the 268 malonylated proteins in the mouse 
liver with two previously acquired malonylation data sets in 
the HeLa cell line by manually searching for protein homologs. 
Our data contained 13 malonylated proteins out of 17 proteins 
identified by Zhao et al. and 54 malonylated protein candi- 
dates out of the 375 proteins identified by Li et al. (supple- 
mental Table S3) (6, 11). Thus, the overlap observed between 
these data sets confirmed the validity of our experimental 
procedures and results. We then compared our data with a 
large-scale succinylation data set, and 135 out of 433 Uni- 
prot-reviewed proteins were shown to overlap between the 
two data sets (supplemental Table S3) (22). These results 
suggested that malonylation and succinylation might occur at 
identical proteins. 

Our data also revealed that some proteins were highly 
malonylated. For example, 29 malonylated peptides (19 
unique malonylated lysine sites) were identified in the one- 
carbon metabolic enzyme 10-FTHFDH. The proteins with 
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more than five malonylated peptides are summarized in sup- 
plemental Table S4. 

Biological Process Analysis of Malonylated Proteins—To 
evaluate the biological relevance of lysine malonylation, the 
PANTHER classification system was used to sort the malo- 
nylated proteins according to their associated biological pro- 
cesses (23). The results of this analysis showed that 194 out 
of 268 malonylated proteins (72%) were classified as involved 
in a metabolic process (Fig. 4D and supplemental Table S5). 
The analysis of Gene Ontology (GO) terms by the Functional 
Annotation Tool of DAVID also revealed that most of the 
enriched GO terms were involved in catabolic or metabolic 
processes (supplemental Table S5) (24). 

Protein Interaction Network Analysis of Malonylated Pro- 
teins—To better understand the function of lysine malonyla- 
tion, the malonylated proteins were subjected to a protein- 
protein interaction (PPI) network analysis using the STRING 
database (25). A PPI network containing 239 nodes and 1517 
edges were constructed. To characterize protein complexes 
among the malonylated proteins, the PPI networks was ana- 
lyzed for highly connected nodes by MCODE (26). Five highly 
connected clusters have been identified including the gluta- 
thione metabolic, oxoacid metabolic, organic acid metabolic, 
glucose catabolic, and translation cluster (Fig. 4E). The PPI 
network analysis indicated that lysine malonylation occurs in 
various metabolic complexes. 

KEGG Pathway Analysis of Malonylated Proteins—To fur- 
ther characterize the role of lysine malonylation, KEGG path- 
way analysis was performed for the malonylated proteins. The 
results showed that malonylated proteins were enriched in the 
glucose, fatty acid, and amino acid metabolic pathways (sup- 
plemental Table S5 and Fig. 5A). Disorders of glucose and 
fatty acid metabolism are two major features of type 2 diabe- 
tes. The majority of enzymes involved in glucose and fatty 
acid metabolism were identified to be malonylated in livers of 
db/db mice (Fig. 5B and supplemental Table S2). Indeed, 
many of the enzymes identified showed elevated levels of 
lysine malonylation (supplemental Table S6). Dysfunction of 
glucose metabolic enzymes, including glucokinase (27, 28), 
glucose-6-phosphatase (29, 30), phosphofructokinase (31), 
and glycogen phosphorylase (32), is known to be involved in 
type 2 diabetes. Although the role of lysine malonylation in 
dysfunction of metabolic enzymes remains unknown, the high 
levels of lysine malonylation on metabolic enzymes in db/db 
livers indicates an association between lysine malonylation 
and metabolic disorders of type 2 diabetes. 

Verification of Lysine Malonylation by Immunoprecipitation 
Coupled with Western blot Analysis—To verify the LC-MS/MS 
results, we selected four proteins (10-FTHFDH, G6PI, FBP1, 
and LDHA) with commercially available antibodies for immu- 
noprecipitation analysis. After immunoprecipitating an equal 
amount of proteins from wt and db/db samples, lysine malo- 
nylation was detected by Western blot. Consistent with the 
LC-MS/MS results, levels of lysine malonylation of 10- 
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Fic. 6. Verification of increased lysine malonylation in db/db 
mice. Four proteins that showed increased lysine malonylation by 
LC-MS/MS method and with commercially available antibody avail- 
able for immunoprecipitation were selected. Tissue lysates were im- 
munoprecipitated by the selected antibodies and protein-A agarose 
beads. For GSTs, tissue lysates were directly immunoprecipitated by 
GSH-beads. Efficiency of immunoprecipitation and malonylation sta- 
tus of target proteins were evaluated by Western blot analysis. A, 
10-FTHFDH (10-formyltetrahydrofolate dehydrogenase). B, Gé6PI 
(Glucose-6-phosphate isomerase). C, FBP1 (Fructose-1,6-bisphos- 
phatase 1). D, LDHA (L-lactate dehydrogenase A chain). E, GSTs 
(Glutathione S-transferases). F-H, Diagram of reactions involving the 
indicate proteins. 10-CHO-FH4 (10-formyltetrahydrofolate), FH4 
(tetrahydrofolate). 


FTHFDH, G6PI, FBP1, and LDHA proteins were elevated in 
db/db livers relative to wt (Fig. 6A-6D). 10-FTHFDH is known 
to catalyze the transformation of 10-CHO-FH4 to FH4 (Fig. 
6F), whereas G6PI, FBP1, and LDHA are important enzymes 
in glucose metabolism (Fig. 6G). In addition, GSTs, which 
catalyze the conjugation of glutathione (GSH) with electro- 
philic compounds, were affinity-purified using GSH-beads 
(Fig. 6H). Elevated lysine malonylation of GSTs in db/db mice 
was detected by Western blot analysis (Fig. 6E). Overall, these 
results verified the elevated levels of lysine malonylation ob- 
served for some of the proteins identified by LC-MS/MS. 
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Fic. 7. Detection of ALDOB activity 
after malonylation of ALDOB and ly- 
sine 107 mutated ALDOB (K107A). 
Change of ALDOB A, and K107A B, ac- 
tivity after in vitro malonylation reaction. 
ALODB and K107A were purified from 
HEK298T cells and incubate at 30°C 50 
with NADH. Activity were determined by 
measuring decrease of absorbance at 
340 nm. p value were calculated based C 
on triplicate measurements. (*: p value < 
0.01, t test). C, Structure analysis of 
mammalian ALDOB protein. The cata- 
lytic center contains several key residues 
including lysine 107 (red). The peptides 
with lysine malonylation (red color K) 
identified by MS/MS are indicated at the 
bottom. 
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Impact of Lysine Malonylation on Enzymatic Activity of 
ALDOB (Fructose bisphosphate aldolase B)—To validate the 
biological role of lysine malonylation, we searched the PDB 
database and found that malonylation can occur on key res- 
idues of several enzymes, including K107 on ALDOB, K124 
on ACAT1 (Acetoacetyl-CoA thiolase), K46 and K273 on 
HMGCS2 (Hydroxymethylglutaryl-CoA synthase, mitochon- 
drial), and K279 and K345 on PBE (Peroxisomal bifunctional 
enzyme) (Fig. 7C and supplemental Fig. S4). To evaluate the 
influence of malonylation on their enzymatic activity, we se- 
lected ALDOB for further analysis. We purified wild type 
ALDOB and K107A-mutated ALDOB proteins and determined 
their enzymatic activity spectrophotometrically. We found that 
enzymatic activity of K107A-mutated ALDOB decreased 30% 
comparing with ALDOB (supplemental Fig. S5). The results 
reveals the importance of K107 for ALDOB activity. We next 
analyzed the influence of K107 malonylation on ALDOB ac- 
tivity. Because of lack of a site-specific malonyltransferase, 
we were not able to obtain K107-specific malonylated 
ALDOB. Alternatively, we nonspecifically malonylated purified 
ALDOB proteins by an in vitro reaction. We found that purified 
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ALDOB was extensively malonylated in vitro, and that enzy- 
matic activity of malonylated ALDOB decreased by 20% 
relative to nonmalonylated ALDOB (Fig. 7A). However, 
enzymatic activity of malonylated K107A-mutated ALDOB de- 
creased little comparing with K107A-mutated ALDOB (Fig. 
7B). Altogether, the body of evidence obtained in our studies 
indicates that K107 malonylation plays an important role in the 
regulation of ALDOB activity. 


DISCUSSION 


In our study, we analyzed acylation patterns in wt and type 2 
diabetic mice. We identified an elevation of protein lysine ma- 
lonylation levels in liver tissues of db/db and ob/ob type 2 
diabetic animal models, and identified large numbers of malo- 
nylated proteins in liver tissues of db/db mice. Bioinformatics 
analysis of the malonylated proteins indicated a role in glucose 
and fatty acid metabolism as their predominant function. In 
addition to liver, we also observed the status of lysine malony- 
lation in other tissues and found that elevation of lysine malo- 
nylation was liver-specific (supplemental Fig. S2). In addition to 
db/db and ob/ob mice, we measured lysine malonylation in 
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other type 2 diabetic animal models including KK mice, ob/ob 
rats, and spontaneously type 2 diabetic monkeys and found 
elevation of lysine malonylation to various degrees (supplemen- 
tal Fig. S3). Thus, we hypothesize that liver-specific elevation of 
lysine malonylation might be common in type 2 diabetes. 

Elevation of lysine malonylation is consistent with previous 
studies reporting higher malonyl-CoA levels in type 2 diabetes 
(12, 13). Because lysine malonylation requires malonyl-CoA, 
the biological function of malonyl-CoA might indirectly reflect 
the biological role of lysine malonylation. Malonyl-trans- 
ferase(s) and demalonylase(s) are also required for lysine ma- 
lonylation. SIRT5 has been identified as a demalonylase (6). 
However, knockdown of SIRT5 only raises lysine malonylation 
moderately, whereas it results in a marked increase in overall 
succinylation levels, suggesting the existence of other, as yet 
unidentified demalonylases (33). In addition, SIRT5 does not 
appear to be associated with type 2 diabetes and SIRT5 
deficiency fails to lead to any phenotypes resembling type 2 
diabetes (33). Future work in this field would have to include 
the identification of the unknown malonyl-transferase(s) and 
alternative demalonylase(s), together with the evaluation of 
their role in type 2 diabetes. 

Similar to lysine malonylation, lysine acetylation has been 
found in almost every enzyme involved in glucose and fatty 
acid metabolism (34). Acetylation of metabolic enzymes is 
known to have a great impact on their catalytic functions (34). 
Acetyl-CoA and malonyl-CoA are closely related and inter- 
convertible, and therefore we propose that lysine malonyla- 
tion also functions in regulating metabolism. In fact, malony- 
lation can occur on key residues of metabolic enzymes (Fig. 
7A and supplemental Fig. S4). Malonylation of ALDOB in vitro 
has been shown to decrease its activity (Fig. 7B). Malonylation 
of other proteins probably affect their functions. The following 
focus in the field would be on studying malonylation of several 
important type 2 diabetic associated proteins. 

Analysis of identified malonylated proteins revealed the role 
of malonylation in metabolic processes. Our malonylation 
data mirror earlier proteomic results for acetylation and suc- 
cinylation, which demonstrated an enrichment of acetylation 
and succinylation in metabolic enzymes (22, 34). Metabolic 
intermediates serve as group donors of protein acylations and 
protein acylations in turn regulate the production of metabolic 
intermediates. This might be a general feedback mechanism 
used in cells to control metabolism in a highly efficient man- 
ner. Given their potential role in metabolic regulation, acyla- 
tions might contribute to different pathological processes. 
Although malonylation might dominate in metabolic disorders 
such as type 2 diabetes, other acylations might be more 
important in other metabolic diseases such as coronary heart 
disease, hypertension, and hyperlipidemia. Recently, eleva- 
tion of lysine acetylation has been detected in the kidneys of 
streptozotocin (STZ)-induced diabetic rats (35). In addition, 
lysine glutarylation has been associated with the disease of 
glutaric academia (36). Altogether, our work identified a large 


number of malonylated proteins, which are potential impor- 
tant in regulating metabolism in type 2 diabetes. 
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